U(l) flux tube profiles from Hamiltonian lattice gauge theory using a 
random walk ground-state projector 
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We use a self-guided random walk to solve the ground-state problem of Hamiltonian U(l) pure gauge theory in 
2+1 dimensions in the string sector. By making use of the electric-field representation, we argue that the spatial 
distribution of the electric field can be more easily measured than in ordinary Monte Carlo simulations. 



1. INTRODUCTION 

Recently, some focus has been given to the ex- 
traction of flux tube profiles from lattice gauge 
theory [1-3]. The usual approach involves cal- 
culating plaquette-plaquette correlators in Monte 
Carlo simulations from which the field distribu- 
tion can be derived. Due to the high level of 
statistical fluctuations, this requires very large 
statistics and has only recently become feasible. 

However, flux tubes and strings arise quite nat- 
urally from the strong-coupling expansion. We 
propose to use the Hamiltonian formulation of 
lattice gauge theory to calculate field-strength 
distributions. The Hamiltonian formulation has 
the advantage of not only having a simple strong- 
coupling limit but can be formulated in the 
electric-field basis, making the field strength a 
fundamental observable. Also, the Hamiltonian 
theory separates in different sectors according to 
string length, reducing the difficulty of extracting 
one particular sector. 

2. METHOD 

2.1. Electric- field representation 

We attempt to solve the ground-state problem 
of the Kogut-Susskind Hamiltonian 



H = 2^ 3 

(ij) 
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Uj is the link operator (electric field) associated 
to thelink ij, U p the plaquette operator (mag- 
netic field) of plaquette p, A ~ 1/g 2 the coupling 
constant. 

In the electric-field representation, the basis of 
the Hilbert space is chosen such that U is diag- 
onal. The basis states of the representation are 
labeled by lattice configurations where each link 
carries a whole number m of electric flux units. 
The plaquette operator (magnetic field) acts as a 
simple raising operator: 
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where n + p represents the basis state n with the 
flux around plaquette p raised by one unit. 
The Hamiltonian reads in this representation: 
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Invariance under gauge transformations implies 
that the Hamiltonian preserves Gauss' law at 
each site, as the total flux at any site is not 
changed by the application of the plaquette oper- 
ator Up. The Hilbert space separates into sectors 
of different external charges that act as sources of 
the electric field. 

In the strong-coupling limit, A — ► and 
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0, and the Hamiltonian is diagonal. In 



each sector, the ground state is the shortest string 
linking the charges. Thus, string states are very 
natural in this model, and numerical simulation 
can center on how they spread out and acquire a 
finite thickness, and how they possibly vanish at 
a phase transition. 

2.2. Random walks 

We employ an ensemble ground-state projector 
using a guided random walk [4,5] to calculate the 
ground-state wave function. This method makes 
use of the similarity between the imaginary-time 
Schrodinger equation 



dt 
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that governs the evolution of a state \ip t ) towards 
the ground-state projector 
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and an ensemble random walk with transition 
probabilities P(n',n) and survival factor V(n): 

A+At(n') = V(n')J2P(n',n)Mn) • 



In an ensemble random walk, members not only 
undergo transitions according to the P(n',n) but 
are also deleted or replicated according to the sur- 
vival factor V(n). 

For a lattice gauge Hamiltonian, ensemble 
members are lattice configurations n. The wave 
function ipt(n) is approximated by the distribu- 
tion of ensemble members in this space. It is 
driven towards the ground-state wave function 
by means of appropriate transition probabilities 
and survival factor P(n', n) and V(n) that are ob- 
tained by equating the discretized imaginary-time 
Schrodinger and the random-walk equation. 

The probability interpretation of P(n',n) and 
V(n) restricts the applicability of the random- 
walk method to certain Hamiltonians. In partic- 
ular, while U{\) lattice gauge theory satisfies this 
restriction, nonabelian theories fail as they would 
require negative probabilities or, equivalently, do 
not ascertain that the ground-state wave function 
is positive. 



2.3. Guidance 

In a many-dimensional problem, the random 
walk easily strays into improbable regions of the 
basis space. Is it thus preferable to offer some 
guidance in the choice of "good" transitions. This 
is done by multiplying the Hamiltonian with a 
guidance wave function $(n): 

H = $" 1 ,ff$ 

(n'\H\n) = —}-^(n'\H\n)$(n) 
$(n') 

The guided Hamiltonian has the same spectrum 
as the original Hamiltonian; its wave functions 
are multiplied ("masked") by the guidance wave 
function. The guidance basically makes moves 
towards the guidance function more probable. 

As the guidance changes the eigenfunctions, 
the extraction of observables requires either to use 
a guidance wave function that is sufficiently close 
to the actual wave function, or to use a split-time 
measurement that removes the guidance. 

We chose a link-separable guidance func- 
tion based on a modified Bessel function and 
parametrized by the center and the width of the 
distribution. As the ground state in the string 
sector is spatially inhomogeneous, different pa- 
rameters can be assigned to each link. We intro- 
duce a form of self-guidance by constantly adapt- 
ing the parametrization to the measured observ- 
ables such that the ensemble described by the 
guidance wave function yields the same observ- 
ables as the actual ground-state ensemble. In 
this sense, the method is variational: it evolves 
the parametrization of the guidance wave func- 
tion towards an optimal parametrization. 

The variational character of this method re- 
duces artifacts from the guidance function as the 
result should be at least as good as one obtained 
from a purely variational method, which would 
be sufficient for measuring the field distribution. 
Is also hints at how the method can be extended 
to systems that are not easily represented as a 
random walk, such as nonabelian theories. 

3. CALCULATIONS AND RESULTS 

We are testing our method in 2 + 1 dimensions 
on U(l) lattice gauge theory. Typically, each 



ensemble has 50 members that are evolved over 
200-300 steps. The ensemble size is held stable 
by a suitable energy offset. While the genera- 
tion process is fast, high statistics (typical 10 6 to 
10 7 configurations) is needed as the wave func- 
tion acquires a large width outside the strong- 
coupling region. The self-guidance introduces an- 
other source of error, as one has to wait until the 
parametrization stabilizes. 
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Figure 1. Electric field strength (in flux units) 
in a cut along the string axis (String of 9 on a 
18 x 11 lattice at A = 1.5). Errors are statistical. 
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Figure 2. Electric field strength in a cut across 
the string axis. 

Fig. 1 shows the distribution of electric field 
strengths in a cut along the string axis for a sam- 
ple string. The deviation from strong coupling is 
manifested in the drop of the field at the center of 



the string accompanied by a spreading out shown 
in the transverse cut in Fig. 2. 

While electric field strengths are easily acces- 
sible in this approach, it is somewhat harder to 
estimate the string tension. The vacuum energies 
dominate the ground-state energies by two orders 
of magnitudes such that calculating the energy 
difference is impracticable. However, we see some 
indication that the classical energy density (E) 2 
may serve as a substitute observable. 

4. SUMMARY 

We presented the application of the guided 
random walk ensemble ground-state projection 
method to the calculation of field distributions 
in flux tubes. This methods targets directly the 
string sector of the theory and thus provides a 
more efficient way to calculate field strengths. 
In particular, they do not appear as plaquette- 
plaquette correlations but as primary representa- 
tion of the wave function. 

In the case of 2 + 1-dimensional (7(1) gauge 
theory, one does not expect a phase transition to 
a nonconfined phase. For further investigations, 
we plan to consider the 3 + 1-dimensional case 
in which a deconfinement phase transition will 
occur. 

To apply the method to nonabelian theories, 
an appropriate scheme is worked at to use the 
random walk as a variational projector which will 
at least return a best-fit positive wave function. 

Finally, inclusion of fermions is straightforward 
in the Hamiltonian formulation both as an ap- 
proximation as well as exactly. 
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